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An important input into reaction tiieory is the density of states or tlie level density. Spectral 
distribution theory (also known as nuclear statistical spectroscopy) characterizes the secular behavior 
of the density of states through moments of the Hamiltonian. One particular approach is to partition 
the model space into subspaces and find the moments in those subspaces; a convenient choice of 
subspaces are spherical shell-model configurations. We revisit these configuration moments and 
find, for complete OTuu many-body spaces, the following behaviors: (a) the configuration width is 
nearly constant for all configurations; (b) the configuration asymmetry or third moment is strongly 
correlated with the configuration centroid; (c) the configuration fourth moment, or excess is linearly 
\ related to the square to the configuration asymmetry. Such universal behavior may allow for more 

^— ■>) . efficient modeling of the density of states in a shell-model framework. 
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I. INTRODUCTION AND MOTIVATION 



An important input into Hauser-Feshbach calculations of statistical capture of neutrons 0,0 is the density of nuclear 
states, or the nuclear level density, as a function of excitation energy. (The phrases "state density" and "level density" 
^ ,. are often used informally as synonyms, although they are not: the state density includes 2 J -I- 1 degeneracies, where 
J is the angular momentum. Because the Hauser-Feshbach formalism uses J-dependent transmission coefhcients, one 
ultimately requires the level density as a function of J. As will become clear below, we discuss both the total state 
QQ ' density as well as J-projected level densities.) 

The density of states is difficult to obtain both experimentally and theoretically and is perhaps the most uncertain 
lO input into Hauser-Feshbach calculations 2]. Therefore there is significant effort, both experimental and theoretical, 
going into determining the density of states. 

The most common theoretical approaches are built around noninteracting Fermi gases, starting from Bethe's 
formula 3] using the single-particle density of states. Of course the residual interaction plays a significant role, 
'-^ '] and single-particle approaches necessarily add corrections to account for the residual interactionj3,|l|. 
^ ■ An alternative is to start from a fully interacting picture such as the configuration-mixing shell model. Interacting 
^ shell-model codes, such as OXBASH 6J, ANTOINE 7], and REDSTICK 8], read in single-particle energies and two- 
r • ■ body matrix elements and compute matrix elements of the many-body Hamiltonian in a truncated, albeit very large, 
many-body basis. All such codes use the Lanczos algorithm to extract the low-lying spectra^- Obtaining the complete 
spectrum is unfeasible for large dimension bases, and unnecessary as well: one only needs the smooth, secular behavior 
^ of the density of states, not the individual states. Instead, one turns to formulations of the interacting shell model 
5^ , that do not involve diagonalization: auxiliary-field path integrals 0, Hit Il2i] , or spectral distribution theory (also 
known as nuclear statistical ST)ectroscoDv [l3l [l4 flSl [l6l| ) . 

Nuclear statistical spectroscopy starts from the moments of the Hamiltonian. For example, it is known that, in a 
finit e sp ace, most state densities of many-body systems with a two-body interaction tend toward a Gaussian shape 
plH Il7| characterized by the first and second moments of the Hamiltonian. In most realistic cases the density has 
small but non-trivial deviations from a Gaussian, so one requires higher moments. It has long been realized that, 
rather than computing higher and higher moments of the total Hamiltonian, one could partition the model space 
into suitable subspaces and compute just a few moments in each subspace. In particular, if one partitions the space 
based upon spherical shell-model configurations, that is, all states of the form (0d5/2)'*(lsi/2)^(0d3/2)^, etc., then it 
is possible to derive expressions for the configuration moments directly in terms of the single-particle energies and 
two-body matrix elements, without constructing any many-body matrix elements (see, however, our discussion in 
depth below). 

Although such an approach has obvious merit, it is far from trivial. For a realistic application, one must include 
thousands of configurations. For example, ^"^Fe in a Is0c?-lp0/ space allowing at most two particles to be excited out 
of the sc?-shell into the pf consists of 147,060 configurations (which in turn represent 9.5 x 10^^ levels). Furthermore, 
higher moments take more and more computational effort. The configuration centroids and widths (first and second 
moments, respectively) take CPU-minutes, while configuration asymmetries (third moments) take CPU-hours or days 
and configuration 4th moments even longer; in our example of ^^Fe it takes about 43 minutes on a single 2.4 GHz 
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Pentium to compute all the second moments, but 9 days to compute the third moments. 

Therefore one important question is: how many moments must be reliably calculated? A related question is: what 
is an appropriate model for the secular behavior? Most work has been based on modified Gaussians, in the Gram- 
Charlier or Edgeware distributions. Kota, Potbhare, and ShenoyT^ made a study of model functions and concluded 
the Cornish-Fisher distribution as the best approximation; they also made a study of the behavior of the configuration 
moments, and this paper can be considered an further investigation of that early work. Zukerflfljl prop osed binomial 
distributions as a way to naturally include asymmetries. Recently Horoi, Kaiser, and Zelevinsky |2(J| have modeled 
the level density as a sum of Gaussians, based only upon first and second moments, and ignoring higher moments 
such as the asymmetry; they argue that the asymmetries tend to cancel. 

Given this background, we have investigated in more detail the configuration moments and their distributions. 
Before sharpening our questions further, we lay out more carefully the formalism of spectral distribution theory and 
configuration moments. 



II. FORMALISM AND MORE DETAILED MOTIVATION 



We begin with the moments of a many-body Hamiltonian, taking some care to define our formalism; this is important 
because we will investigate if the behavior of moments depends upon, for example, angular momentum. None of the 
formalism in this section is original; a thorough reference to statistical spectroscopy is Ref. (16j , although some of our 
notation is different. 

We always work in a finite model space wherein the number of particles is fixed. If in Ai we represent the 
Hamiltonian as a matrix H, then all the moments can be written in terms of traces; for some matrix O the trace is 



trO= 5](z|0|z) (1) 



ieM 

The total dimension of the space is Z? = tr 1, and then the average is 

(0) = ltrO. (2) 

The first moment, or the centroid, of the Hamiltonian is then 

^ = (H); (3) 
all other moments are central moments, computed relative to the centroid: 

= ((H - ^)") . (4) 
The moments can also be computed directly from the density of states 

p(£;) =tr^(£;-H), (5) 
E = D-'^ j dEEp{E), (6) 

= D-^ J dE{E- EYp{E). (7) 
The width cr is given by \/JiJ^, and one typically scales the higher moments by the width: 

m(«) = ^. (8) 

In addition to the centroid and the width, the next two moments have special names. The scaled third moment 
to'^' is the asymmetry, or the skewness, while m^^^ — 3 is the excess, so called because for a Gaussian m^^^ — 3 and 
hence a Gaussian has zero excess. 

In all of the above, one has to be careful to clearly state the model space Ai. If, for example, M contains only 
states of a fixed Jz, then one has the level density. On the other hand, if Ai contains all values, so that one has 
all 2 J -I- 1 degeneracies, one has the state density. 

Next we consider projection operators and subspaces. For this paper we use subspaces based upons spherical 
single-particle configurations, that is, sets of states of the form, e.g., (0c?5/2)^(0d3/2)^, (0^5/2)^(0^3/2)^(151/2)^, etc.. 
Angular momentum alsos play an important role, so we adopt nomenclature as follows: 
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A J-summed configuration means a sum over all J and all Jz in the configuration; J-projected means a sum just 
over states of fixed J in a configuration. We can relate the J-summed density, which is really the state density, to a 
sum over J-projected (level) densities: 

p{E)^Y.^2J+l)pj{E). (9) 
J 

One can in principle also project out states of good isospin. For simplicity we consider only pn-configurations, 
that is, configurations that have the proton and neutron occupations fixed separately, i.e., (0^5/2)^(0^3/2)^ is distin- 
guished from (0^5/2)^(0^3/2)^ (which, for an isospin-conserving interaction, will have the same moments) and from 
(0^5/2)^(0^5/2)^,(0^3/2)^(0^3/2);', (which, for an interaction that does not mix isospin but does depend on isospin, 
will generally have different moments). We do not believe that our results will differ dramatically for projection onto 
isospin. 

Finally we turn to the key question of subspaces. We use a, /3, 7, . . . to label subspaces. Let 

Po.^Y.\')^\ (10) 

be the projection operator for the a-th subspace. One can introduce partial or configuration density^ 

p^{E)=ivPJ{E-n). (11) 

The J-summed and J-projected configuration densities are defined in the obvious way (and because projection of 
angular momentum commutes with projection onto single-particle configurations, leads to no difficulties), and can be 
related by pc,{E) = Ej(2 J + l)pc...j{E). 

Now we can define configuration moments: the configuration dimension is Da — trPa, the configuration centroid is 

Ea^D-hiP^H, (12) 

while the configuration width Oa is given by 

al^D-hrPaiH-Eaf, (13) 

the configuration asymmetry is 



m 



(3) = ^i?-itrP„(H-^„)3))/a3, (14) 



and the configuration excess is defined similarly. 

Note that in all the above traces, one inserts only one projection operator. This means that, although one is working 
with the "density" in a particular subspace, matrix elements of the Hamiltonian can take one out of the subspace. 
For example, one can further define, distinct from the configuration width, the partial variance (sometimes called, 
misleadingly, the partial width) 



J. 1 

i a/3 = 



D 



73 - Oal3 



-^tr (PaHP^H) - 5^pEl 



(15) 



Then the configuration variance a\ arises from two contributions, an intrinsic piece V aa and a sum over extrinsic 
pieces V^fi- Although some authors make extensive use of partial widths, we will not use them in this paper except 
to understand the origin of correlations in the configuration asymmetries. 



A. Monopole and traceless interactions 



One useful result from spectral distribution theory is that the configuration centrmds depend entirely upon the 
single-particle energies and the monopole- monopole part of the residual interaction(2]|. The monopole interaction 
is attributed to mean-field and saturation properties of the nuclear interaction. One can subtract out the monopole 
interaction exactly, which sets the centroids to zero but leaves the widths unchanged; such a monopole-subtracted 
interaction is referred to as a "traceless" interaction. Third and fourth configuration moments are changed by 
subtracting out the monopole interaction, as we will see below. 

We briefly review calculation of the centroids because they give the nonexpert a better idea of the methods of spectral 
distribution theory. We label spherical orbits such as 0(^5/2, lsi/2, etc, by a,5, c, .... Let Eq be the single-particle 
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energy for orbit a and let \ah,J) be an antisymmeterized, normalized 2-particle state with angular momentum J (it 
turns out that the label is unnecessary). Then let Vj{ah,cd) be the matrix elements of the Hamiltonian between 
|a6, J) and |cd, J), after subtracting off the single-particle energies. The monopole potential is then given by 

iVQ(iVf, +dab) ^ 

where Na = 2ja + 1 is the degeneracy of the orbit a. 

Any spherical configuration a is defined by the occupation of the spherical orbits ni, 712, J^s, • ■ •• Then the (J- 
summed) configuration centroid is 

Ea = ^riflCa + ^^"-a("-6 " 5ab)U{ah). (17) 

a ah 

A traceless interaction is made by subtracting off the monopole potential, AVj {ab, cd) = Vj {ab, cd) — dacSbdU (ah) . One 
uses the traceless interaction to compute the configuration width, although the resulting formula is now fourth-order 
in the Ua- 

One can extend the above formulas to isospin, although we find it convenient to work in a pn formalism where 
protons and neutrons are assigned to distinct orbits. 



III. SOME QUESTIONS 

Now that we have reviewed the required mathematical language we can ask some interesting questions. For example: 
Ql: What is the distribution of the configuration widths? Do they depend upon the dimension of the configuration? 

Upon A? We find: the configuration widths are remarkably constant for a given nuclide, but depend upon A and 

N - Z. 

Q2: What is the distribution of the configuration asymmetries? What range do the configuration asymmetries 
have? To model the partial (configuration) level densities it may be useful to know the range of the configuration 
asymmetries rria , since many numerical models often have restrictions in this respect; for example the binomial 
distribution [13 has difficulty with Im'^'l > 1 p3|. 

Are the asymmetries correlated with, for example, the configuration centroid? This is relevant because Horoi, Kaiser, 

and Zelevinsky 20] speculate that the asymmetries in configuration densities tend to cancel out. This speculation is 

ill-founded, because in fact there is a strong correlation of the configuration asymmetry with the the configuration 

centroid. In particular, configurations low in energy often, but not always, have a significant positive asymmetry 

(skewed towards higher energies) while configurations at high energy have significant negative asymmetry. (Such 

a correlation can be seen in Fig. 1 of Ref. [T^l, although those authors did not comment on it.) This means, at 

any given excitation energy, all of the "nearby" configurations that contribute to the level density have nearly the 

same asymmetry. We note, however, that Horoi et a^|20|| claim satisfactory modeling of the level density without 

asymmetries; while this result requires further investigation, we speculate they may have happened upon fortuitous 

cases where the asymmetry at low energy is small. 

The correlation between the asymmetries and the centroids is intriguing. We suggest that a large asymmetry arises 

(3) 

through connections to distant configurations. One can divide the configuration third moment rria into two parts, 
an intrinsic piece which comes strictly from matrix elements that do not take you out of the subspace, and a part 
that looks like 

Y,Dpr^0{Ep-E^). (18) 

13 

Such terms could explain the correlation between the configuration asymmetries and the configuration centroids. We 
test this idea by setting all the centroids Ea = 0, that is, by considering so-called "traceless" interactions. When we 
do this, we find the configuration asymmetries are dramatically reduced. Furthermore, deviations from this trend we 
ascribe to collectivity. 

Q3: What is the distribution of the excess? In particular we investigate correlations between the excess (fourth 
moment) and asymmetry (third moment). If there were a correlation, it would help in modeling the state density, both 
in the choice of appropriate model functions for the secular behavior, but also save computations: one could compute 
only the asymmetry and estimate the excess. We find a strong correlation, that in fact the excess is linearly related to 
the square of the asymmetry, with high statistical significance. The exact parameters of the linear correlation depend 
weakly upon the interaction and the nuclide in question. 
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FIG. 1: (Color online) Distribution of configuration in *'*Ti as function of the corresponding configurations centroids, in the 
p/-shell with interactions GXPFl ^ and FPD6G 



IV. DETAILS OF CALCULATION METHODS 



We work in a spherical shell-model basis, and the input to our calculations are the same as that of most spherical 
M-scheme codes: a list of single-particle orbits, and a list of single-particle energies and two-body matrix elements. 
In particular we work in the lsi/2-0(i3/2-0d5/2 or sd space, where we use the Brown- Wildenthal universal sd (USD) 
interaction^], and in the Ipur^P^/rOfj^./rOfr/?. or pf shell, where we use the interactions GXPFl [2^, FPD6G 
[2^|2^, and KB3G[2^ |23, Hg- We also used quadrupole-quadrupole {Q ■ Q) interaction and Gaussian-distributed 
random two-body interactions (the two-body random ensemble or TBRE) in both spaces to investigate any dependence 
of our results on the interaction. Finally, we also create "traceless" interactions, by subtracting off the monopole 
potential, Eqn. \i<6\ . so that the centroids = identically. 

In principle, one can compute J-summed, pn-configuration moments of order one through four, starting with single- 
particle energies and two-body matrix elements, without constructing the many-body matrix elements^lS, .14, 16]. We 
have found the formulas for first and second moments work as stated, but however we were unable to satisfactorily 
validate the trace formulas for third and fourth configuration moments. (At this time we are still unsure if the fault 
lies in the published formulas or in our implementation of them.) 

Therefore for this work wc used the REDSTICK shell- model codeQ (which uses input files similar to those for 
OXBASHj^]) to construct the many-body matrix elements and to compute the configuration moments directly. This 
is inefficient, and clearly for large-scale computations one would need to satisfactorily implement trace formulas for 
all required moments. On the other hand, our direct method allows us to exactly compute J-projected moments as 
well, something which cannot be done exactly by trace formulas, although approximate methods do exist 16, 29, 3Qj]. 

For this paper we consider only complete Ohw many-body spaces, due to our current computational limitations. It 
will be interesting in future work to consider multi-shell spaces. 



V. RESULTS 



A. Distribution of configuration widths 



We begin by considering the distribution of configuration widths. It is already known that the configuration 
widths are approximately constant [l^. which we illustrate in Fig^for '''^Ti in the p/-shell (2 valence protons and 2 
valence neutrons), using all 100 J-summed pn-configurations. We consider two interactions, GXPFl (filled circles) 
and FPD6G (open circles), which, unsurprisingly, yield similar results. The spreading of the widths corresponding to 
the FPD6G interaction is very narrow. The GXPFl widths show a slightly greater range. 

While the configuration widths are constant for a given nuclide, they do depend on the nuclide in question. Ta- 
ble shows the average configuration widths (and corresponding standard deviations) for a number of representative 
nuclides with valence particles in the sd-shell with the USD interaction '23']. The widths are largest for "half-filled" 
cases and for TV ~ Z nuclides. The scatter, as represented by the standard deviations, is nearly constant as well. 
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TABLE I: Average configuration widths and corresponding standard deviations for several nuclides in the sci-shell with USD 
[2^ interaction. 
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B. Distribution of configuration asymmetries 



In Fig. 121 we plot the configuration asymmetry against the configuration centroid for two representative examples, 
■^■^Ar in the sc?-shell with the USD interaction (this is similar to Fig. 1 of Ref. for 8 particles in the sd-shell, also 
showing rather strange gill-like structure), and '*'*Ti in the p/-shell using GXPFl. Other nuclides in these spaces show 
the same generic behavior. The key results are, first, the range of the asymmetries, and second, the correlation of 
asymmetries with centroids. Please remember that here and throughout we consider scaled asymmetries and excesses, 
that is, normalized by the relevant configuration width. 

The configuration asymmetries range from -2 to 2, which seems to be a general limit; in other sd- and p/-shell 
nuclei, we found only rare cases of \ma | > 2 for J-summed moments ( J-projected moments have wider ranges as we 
will see in a later section). It is important to keep in mind that these findings are only for complete Ohfl spaces. We 
have not yet considered multi-shell configuration which may have a different range; such spaces are beyond the scope 
of our current methods and must be investigated in future. 

The asymmetries are strongly correlated with the configuration centroids Ea'- the most positive asymmetries are 
found at the most negative centroids and the the most negative asymmetries at the most positive centroids. 

Following the discussion in lllll we recomputed the asymmetries for "traceless" interactions (open squares) in Fig|21 
that is, we subtract off the monopole potential. (In such a case the centroids all automatically go to zero, but we plot 
the asymmetries against the original centroids so one may see how the asymmetries are shifted.) Now most of the 
configuration asymmetries go to nearly zero. This supports the hypothesis that much of the configuration asymmetry 
derives mostly, although not completely, from Eq. (fl^ . 

We considered a number of nuclides and interactions and found consistently the generic behavior illustrated in 
Fig. 121 There are two important variations. First, we refer the reader to Fig. 1 of Ref 18], which for 8 particles in the 
sd shell has a range of asymmetries from 0.5 at low energies to -2 at high emergies. As another example, in Fig. Olwe 
show the distribution of configuration asymmetries for ^''Ne, using the USD interaction. Several of the asymmetries 
for low-lying configuration, in particular the lowest energy configuration, have asymmetries close to zero. While these 
are important exceptions to the examples in Fig. [21 the converse is also true: suggestions that asymmetries are small 
at low energies 0, do not generalize. 

For further investigation, in Fig. Owe again consider the asymmetries for a traceless USD interaction in ^°Ne. Most of 
the asymmetries again go to zero, but the low-energy asymmetries that were near zero now become large and negative. 
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FIG. 2: (Color online) Configuration asymmetries versus centroids for '^'^Ar (top) and **'Ti (bottom). Filled circles are for the 
full interaction; open squares have the monopole part of the interaction removed ("traceless"). 
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FIG. 3: (Color online) Configuration asymmetries versus centroids for ^"Ne. Filled circles are for the full interaction; open 
squares have the monopole part of the interaction removed. 



We postulate that this arises from (probably attractive quadrupole) collectivity in the lowest states, which would lead 
to a large negative asymmetry. To test this idea we plot in Fig. 0] asymmetries for a pure (repulsive) Q ■ Q interaction, 
both with and without ( "traceless" ) the monopole- monopole interaction. In all our cases, including the two examples 
from the so? shell ('^'^Ar, top) and pf shell ('*'*Ti, bottom) we see strong asymmetries, even for the traceless interaction. 
In hindsight this is not surprising for a collective interaction, but supports our interpretation of the behavior of more 
realistic interactions in Figs. |(2J, Q. A realistic interaction has both attractive and repulsive collective pieces, e.g. 
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FIG. 4: (Color online) Configuration asymmetries versus centroids for the collective Q ■ Q interactions for ^^Ar (top) and "^^Ti 
(bottom). Filled circles are for the full interaction; open squares have the monopole part of the interaction removed. 



Q ■ Q and S ■ S, respectively, and in the bulk the coUectivity tends to cancel, leading to traceless asymmetries near 
zero. At low energy, however, the attractive quadrupole collectivity dominates, pushing the asymmetry lower. This 
effect is smaller in the pf shell, where a larger spin-orbit splitting weakens quadrupole collectivity |3lj |. To summarize, 
the behavior of the configuration asymmetries and any correlation with the configuration centroids result from an 
interplay between collective interactions and the monopole-monopole interaction. 



C. Correlation of asymmetry and excess 

We now turn to look at the relationship between the third and the fourth moments. Figure 1 of Ref. 'I8l indirectly 
suggest a correlation but those authors did not comment on it. Fig. [Slplots (ma^'')^ versus m'^'^ for two nuclides in the 
sc?-shell valence space, ^"Ne and ■^^Ar. The observed linear correlation occurs in all other nuclides we considered as 
well, and Table UTI gives the slope and intercept (as well as the corresponding "uncertainties," which are the diagonal 
correlation coefficients) for several sd-shell nuclides, all with the USD interaction|^. Table ITlTl shows slope, intercept, 
and uncertainties for a single nuclide in the p/-shell, ^"^Ti, but for a variety of interactions, including quadrupole- 
quadrupole and a member of the two-body random ensemble (TBRE). There is a weak but statistically significant 
(the differences are several standard deviations) dependence on interaction and nuclide. 

We attribute the relation between m'-'^) and m^^^ to a maximum entropy principle. For convenience we took the 
distribution a modified Breit-Wigner, 
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lines are least-squares linear fit. 



TABLE II: 


Intercept, slope, and respective 


standard deviations (era, o";, 


corresponding 


to linear fits of m4 vs. 


m| plots for 


different nuclides in the sd-shell valence space, with the USD interaction 


El. 






Nuclide 
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^intercept 
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33 Ar 
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0.92 




0.036 


0.033 



and can adjust the parameters to fix the moments. We chose the form (|19|l because it is positive definite on the 
interval (£'rnin, ^-max), and one can compute the moments analytically. (Some alternatives such as Gram-Charlier are 
not positive-definite; for yet other alternatives, such as Cornish-Fisher 18] and binomial [T9| . in practice one must 
integrate numerically to get accurate moments at large asymmetries.) We then computed the entropy, 

S = -[ p{E)\np{E)dE. (20) 

For fixed width and normalization, the maximum entropy describes a straight line in the m^'^^)^-™^'*^ plane-albeit 
with different parameters, a slope of 1.67 and an intercept of m'--^^ — 3 = —0.07. Nonetheless, this simple model makes 
for a reasonable understanding of the correlation between the third and fourth moments. 



D. Distribution of J-projected moments 



Finally, we consider the behavior of J-projected moments. This is relevant because some authors advocate the use 
of symmetric (Gaussian) distributions for J-projected partial densities, eschewing moments higher than the width 



TABLE III: Same as Table ITTI but for different interactions on p/-sheU "''Ti. 
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FIG. 6; (Color online) Comparison of J-projected moments for ■*''Ti for two configurations. 



In Fig.Elwe plot J-projected moments for two representative configurations for "^^Ti using GXPFl. One configura- 
tion (the 100th partition, or PlOO) is nearly symmetric with = —0.12, while the other (P079) has high asymmetry 
(ma = 1.02). The centroids are plotted relative to the corresponding J-summed centroid, for ease of comparison. 

The moments are relatively constant, although with larger scatter; this not surprising because of the smaller 
dimensions, and the reader is cautioned not to overinterpret these graphs. (For partition 100, the total dimension is 
784 levels, with 4 J = levels, 130 J = 6 levels, and 25 J = 12 levels; for partition 79, the total dimension is 1536 
levels, 4 J = levels, 234 J — Q levels, and 23 J = 11 levels.) The most significant trend is the width, which is largest 
for J = and shrinks almost monotonically. This was found in previous investigations for J-projected widths over 
the entire spacels^, 13 been related to the predominance of J = ground states even with random two-body 

interactions '35', '36!| . Similar behavior is found in other configurations in '^'^^i as well as in other nuclides. 

In Fig. 13 we compare asymmetries versus centroids for both J-summed and J-projected moments (again for '*'*Ti). 
The J-projected moments track the J-summed moments. 

Finally in Fig. |Slwc consider the correlation of TO4 against (tos)^ for both J-projected and J-summed moments 
(again for ^^Ti). Both show the same linear relation, although with slightly different slope and intercept. Also, the 
J-projected moments have a greater range. 



VI. SUMMARY 



We have investigated in detail the distribution of and correlations among configuration moments in QTioj sd- and 
light pf-sheW nuclei. We regained earlier results finding: 

• the configuration widths are constant (but depend on the nuclide in question); 

• J-projected widths are largest for smallest J and decrease as J increases; and 

• the configuration asymmetry (third moment) is correlated with the configuration centroid (which, although docu- 
mented before in the literature, had never been commented on). 

We also found some new results, namely: 

• the configuration asymmetries arise mostly, though not entirely, through the monopole part of the interaction; and 

• the configuration excesses are linearly dependent upon the square of the asymmetry, which we attribute to a 
maximum entropy principle. 

These last two are particularly interesting. They suggest that, if one wants to include the third and fourth moments, 
which are computationally intensive, one might be able to take a short cut and approximate the asymmetries from 
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FIG. 7; (Color online) Asymmetries versus centroids for ''''Ti in the p/-shell with interaction GXPFl i24l . We compare 
J-summed asymmetries (filled circles) against J-projected (triangles). 




FIG. 8: (Color online) Excess (mi) plotted against the square of the asymmetry (ma) for *'*Ti in the p/-shell with interaction 
GXPFl M. We compare J-summed asymmetries (filled circles) against J-projected (triangles). Straight lines are least-squares 
linear fit. 



Eqn. (|18() and the fourth moments from the Hnear relationship with the square of the asymmetry. Neither "universal- 
ity" is absolute; intrinsic asymmetry (collectivity) at low energy must be considered to get correct asymmetries, and 
the slope and intercept of the m^'^^-m'-'*) relation depends weakly upon the interaction and nuclide. In the near future 
we plan to investigate in detail how much error is induced through such approximations and indeed if using the third 
and fourth moments do or do not improve descriptions of the density of states. 

There is one final consequence of our work: if one is indeed to exploit knowledge of the asymmetry, any "model 
function" must be able to handle a range of at least —2 < rria^ < 2. For example, binomial distributions |0| fail to 
work over this entire range. We will address this issue in future work. 
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